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Abstract 

This talk describes the resummation approach in (Fractional) Analytic Perturbation 
Theory (FAPT) in QCD. First, we make a short historical review of the (F)APT approach 
and then shortly describe the global scheme of FAPT which allows one to take into account 
heavy-quark thresholds. After that we show how it is possible to resum a non-power series 
in (F)APT both in the one- and two-loop approximations. As an application we suggest 
our analysis of the Higgs boson decay H° — > bb, important for the LHC program, and of the 
vector-current Adlcr function. 
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Cl; 1 Analytic Perturbation Theory in QCD 
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First, I'd like to remind the history of Analytic Perturbation Theory (APT). Strictly speak- 
ing, APT was initiated by the paper of N. N. Bogolyuov et al. in 1959 [lj, where ghost-free 
effective coupling for QED was constructed using the dispersion-relations method. Then, in 
1982, Radyushkin [2], and Krasnikov and Pivovarov [3], using the same dispersion-relations 
technique, suggested regular (for s > A 2 ) QCD running coupling in the Minkowski region, 
a s (s) = a s (s)/Pf, namely, a s [L] = ir^ 1 arctan(-7r/L) with L = ln(s/A 2 ) and (3f = bo(N f) / (4tt) = 
(11 - 2iV//3)/(47r)Q (This normalized coupling a s [L] in standard QCD PT is just l/L in the 
one-loop approximation.) 

Thirteen years later, Dokshitzer et al. constructed an IR finite af ff (Q 2 ) in the Euclidean 
domain using a renormalon-based approach with dispersion relations [5J. After that, Jones and 
Solovtsov [3 [6] discovered the coupling which appears to be finite for all s and coincides with 
the Radyushkin's one for s > A 2 , namely 

2ti[L] = / da = — arccos = . (la) 

Js ° n V7T 2 + L 2 

Just at the same time, Ball, Beneke and Braun reproduced in the pQCD-based renormalon 
approach a s [L] = tt^ 1 arctan(7r/L) in the low-energy Minkowski region 0[H], and Shirkov and 
Solovtsov [9], using the same dispersion-relations approach of [1], discovered the ghost-free 
coupling, 

^'T^-i-^i- (lb) 



* E-mail: bakulev@theor.jinr.ru 

1 We use the notations f(Q 2 ) and f[L] in order to specify the arguments we mean — squared momentum Q 2 
or its logarithm L = ln(s/A 2 ) or ln(Q 2 /A 2 ), that is f[L] = /(A 2 • e L ) and A 2 is usually referred to the N f = 3 
region. 
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in the Euclidean region with L = ln(Q 2 /A 2 ). But Shirkov-Solovtsov approach, named APT, 
appears to be more powerful: in the Euclidean domain, — q 2 = Q 2 , L = lnQ 2 /A 2 , it generates 
the set of images for the effective coupling and its n-th powers, {-4n[£]} ne pji whereas in the 
Minkowski domain, q 2 = s, L = lns/A 2 , it generates another set, {2ln[-^]}„ e N (see also in [10J). 
APT is based on the Renormalization Group (RG) and causality that guarantees standard 
perturbative UV asymptotics and spectral properties. The power series d m a™[L\ transforms 
into a non-power series ^ m d m 2t m [L] or ^ m d ra i m [I] in APT. The global version of APT, 
which takes into account the heavy-quark thresholds, has been elaborated in [llj . 

By the analytization in APT for an observable f(Q 2 ) we mean the "Kallen-Lehmann" 
representation 

with Pf(cr) = — Im [/(—a)], which transforms to different prescriptions for couplings in the 

7T 

Euclidean and Minkowski domains 

ML] = A E K]= H-^-da; (2b) 
Jo ° + v 

ML\ = A M K]= r^-da. (2c) 

Js & 



In the one-loop approximation, we have then pi(a) = L 2 + ir 2 and Eqs. ([I]), whereas 
analytic images of the higher powers (n > 2, n G N) are [2j ?] 



ML}\ 1 ( d \ n_1 / A\[L] 



ML] J (n-1)! \ dLJ \%i[L 



(3) 



At first glance, the APT is a complete theory supplying tools to produce an analytic answer 
for any perturbative series in QCD. But in 2001 Karanikas and Stefanis reminded the need 
for more tools to produce analytic expressions for hadronic observables, calculated perturba- 
tively [12] , Indeed, in the standard QCD PT one has also: 

(i) the factorization procedure in QCD that gives rise to the appearance of logarithmic factors 
of the type: a%[L] L; 

(ii) the RG evolution that generates evolution factors of the type B{Q 2 ) = [Z(Q 2 )/Z(p 2 )] 
B(p 2 ), which reduce in the one-loop approximation to Z(Q 2 ) ~ a v s [L] with v = 70/(260) 
being a fractional number (here the subscript o means that the corresponding quantity is 
calculated in the one- loop approximation). 

They suggested the principle of analytization "as a whole" in the Q 2 plane for hadronic am- 
plitudes, obtained in QCD PT, thus generalizing the analytic approach suggested in [13] I A 
review of both approaches can be found in [TJ], whereas an extensive discussion of the APT 
application to the pion form factor analysis is contained in |15| . 

All that means that in order to generalize APT in the "analytization as a whole" direction 
one needs to construct analytic images of new functions: a", a v s L m , .... This task has been 
performed in the framework of the so-called FAPT, suggested in [16^ |17j . Now we briefly 



2 More precisely, they proposed the analytization recipe for terms like J^dxj^dya B {Q 2 xy) f(x)f(y), which 
can be treated as an effective account for the logarithmic terms in the next-to-leading-order approximation of 
perturbative QCD. 
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describe this approach. In the one-loop approximation, using the recursive relation ([3]), we can 
obtain explicit expressions for ^[L] and 21. [L]: 

a rn 1 F{e- L ,l-u) 

A ^ = T»~ W) ; (4a) 



sin 



(^-l)arccos(^A= 



*» W = 77—, . ^:Tw/ J • (4b) 



vr(^-l)(^ 2 +L 2 ) ( ^ 1)/2 

Here F(z, v) is the reduced Lerch transcendental function, which is an analytic function in v. 
This function has very interesting properties, which were discussed extensively in the previous 
papers QSJ EJ Q21 02] . 

The construction of FAPT with a fixed number of quark flavors, Nf, is a two-step procedure: 
we start with the perturbative result [a s (Q 2 )] , generate the spectral density using p v {cr) = 

— Im \a s (— a)] " , and then obtain analytic couplings *4..[L] and 21. [L] via the relations like Eqs. 

7T 

([2]). Here Nf is fixed and factorized out. We can proceed in the same manner with A^-dependent 
quantities: [a s (Q 2 ; N f )] u =► p v {a;N f ) = p v [L a ;N f ] = p v {a)/P v f => A u [L;N f ] and %,[L;N f \ 
- here Nf is fixed, but not factorized out. The global version of FAPT [18], which takes 
into account heavy-quark thresholds, is constructed along the same lines but starting from the 
global perturbative coupling [a| lob (Q 2 )] , being a continuous function of Q 2 due to adopting 
different values of QCD scales A/, which correspond to different values of Nf. We illustrate 
here the case of only one heavy-quark threshold at s = m 2 , corresponding to the transition 
Nf = 3—7- Nf = 4. Then we obtain the discontinuous spectral density 

pf » = 6{L a < L A ) p n [L a ; 3] +9 (L 4 < L a ) p n [L a + A 4 ; 4] (5) 



with L a = In (cr/Af), Lf = In ym 2 f/ls^j and Xf = In \^A 2 /A 2 J for / = 4, which is expressed in 

terms of fixed-flavor spectral densities with 3 and 4 flavors, ~p n [L; 3] and ~p n [L + A4; 4]. However, 
it generates the continuous Minkowski coupling 

2lf b [L] = 6{L<U) (%[L;3] + A 43 %) +9(L 4 <L) 2T l/ [L + A 4 ;4] (6a) 



with A432l„ = Stj, [L4 + A4; 4] — 21^ [L4; 3] and the analytic Euclidean coupling ^ lob [L] 



AriL] =A[L + A 4 ; 4] + / ^L^ -jAL a + \,-A] dL ^ 

1 + e CT 



The same strategy one needs to use when working in a higher-loop approximation — then, 
analytic formulas even for the fixed- Nf case are unavailable (the spectral density is expressed 
in terms of the imaginary part of the Lambert function) and one is forced to use integral 
representations of the type ([2]) to calculate analytic couplings and other quantities of interest. 
One more complication here is the more complex evolution law: for example, in the two-loop 
approximation the evolution factor B(Q 2 ) = [Z (Q 2 ) / Z (p 2 )] has Z{Q 2 ) ~ a£°[L](l + c 1 a s [L]) V1 , 
with exponents vq and v\ being known numerical coefficients — for more details see in |18| . 



2 Resummation in (F)APT 

We consider now the perturbative expansion of a typical physical quantity, like the Adler func- 
tion and the ratio R, in the one-loop APT. Due to the limited space of our presentation, we 
provide all formulas only for quantities in the Minkowski region: 

00 

TZ[L] = J2dn% l [L]- (7a) 

n=l 
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We suggest that there exist a generating function P(t) for the coefficients d n = d n jd\: 



OO |>00 

d n = I P(t)t n ~ l dt with / P(t)dt = l. 
lo Jo 



(7b) 



To shorten our formulae, we use for the integral / f(t)P(t)dt the following notation: ((f(t)))p(t) ■ 

o 

Then, the coefficients d n = d\ {{t' 11 ^ 1 )) pm and the series ([Taj) can be represented as follows (we 
put for shortness d\ = 1): 



K[L} = J2((t n - 1 ))p(t)%n[L}. 



(7c) 



n=l 



2.1 One-loop case 

Due to the recurrence relations 
sum in ATcl). viz., 



as has been shown in j20J, we have the exact result for the 



K[L] = {{% l [L-t\)) P{t) 



(8) 



The integral over the variable t has a rigorous meaning, ensured by the finiteness of the coupling 
2ti[t] < 1 and the fast fall-off of the generating function P(t). This result resembles the result 
of resummation of the large (/3q a s ) n ~ 1 terms due to Neubert [21] . but here we have no problem 
with renormalon poles, i. e., in our case the integral ((2ti|X — t])) pm is nicely convergent. Note 
that the same type of approach in the one-loop APT has been independently invented two years 
lateJl by Cvetic and Valenzuela [22J within the framework of the so-called "skeleton expansion" . 

We have first accomplished the generalization of result ([S]) to the case of global APT [11] , 
taking into account heavy-quark thresholds [HI [23]. Then, one starts with the series of type 
([7c]) , where 2l n [L] are substituted by their global analogs 2l^ lob [L] (note that due to the different 
normalizations of the global couplings, 2l^ lob [L] ~ 2l n [L]//3y, the coefficients d n should also be 
changed). Then, we have 



K sloh [L] = 0(L<L 4 )((A 4 2li[t] +21] 



L- 



)) P(t) + e(L>u)((% 



L + \ A -— ;4 

P4 ■ 



P(t) , 



(9) 



where A 4 QL v [t] = 2l„ L 4 + A 4 - t/fc 4 -%,L 3 - t//3 3 ; 3 

The second generalization has been obtained for the case of global FAPT. The starting point 
in this case is the generalized recurrence relation 



"77 J~n+u [L] — J- n +i+ u [L] , 



n + v dL 



(10) 



where J-[L] denotes one of the analytic quantities A[L] or 21 [L]. The result of summation of the 
series Yl^Lodt 11-1 )) p(t) ^n+v[L] is the complete analog of Eq. Q with the substitutions 



where d =>■ do% v [L\, 2ti[L - 1] 2li +i ,[L - t], and A 4 2ti[t] => A 4 2li +i/ [t]. All needed formulas 
have also been obtained in parallel for the Euclidean case, see in 



3 Note that the paper of Mikhailov [2D] was put into arXiv in November, 2004, and after a long discussions 
with referees published only in 2007, whereas the paper of Cvetic and Valenzuela [22] was put in arXiv in August, 
2006, and published the same year. It is interesting to remind here the title of Mikhailov's paper: "Generalization 
of BLM procedure and its scales in any order of pQCD: A practical approach" — we see that it was devoted to 
completely different subject and Eq. (|8|> , numerated in this paper as Eq. (2.7), was obtained in passing. 
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2.2 Two-loop case 

In the case of the two-loop running of cg\Q 2 ) = An a^(Q 2 )/b , i-e., when the /3-function has 
the following form (with c\ = &1/&0) 

/3(«(2)) = -«( 2 ) 



[l + ci a {2) ] 



(12) 



we have two types of complication. First, the recurrence relations ([3]) and (jlOp transform into 

1 d 



(13) 



n + z/ cZL 

where .F[L] denotes now one of the analytic quantities „4.( 2 )[L], 2l( 2 )[L], or p^[L]. Second, the 
evolution factors appear to be more complicated as well, namely, 

Z U0 ., Ul [L} = a^[L}{l + c 1 a {2) y i [L], (14) 

with exponents vq and v\ being numerical combinations of anomalous-dimension coefficients 70, 
71, and /3-function coefficients /3q, f3%. 

For these reasons, we need to consider the resummation of the following series: 



s v [L] = ^((^ 1 ))pw^W; 



n=l 

00 



J v ,v 1 



[L] 



^(r- 1 )) P(t) J-, 



[L] 



(15a) 
(15b) 



n=l 



(2) (2) 

In the last line, one has T n +vo,vi [L] = Ul [L] or 25„_j_j, ^ IAI ( or — m the case °f global 

(2) 

FAPT — the spectral density p n + uo vx [L] for these functions) with 



Bv Ml [L] 



A E 
A M 



l {2) 



,1 + ci a {2) ) 



[L] 



(15c) 



being the analytic images of the two-loop evolution factors in the Euclidean and Minkowski 
regions, correspondingly. These problems have been successfully resolved in our recent paper [23] 
and here we show only the results. 

First, we need to define the two-loop evolution "logarithmic time" r(i): 



r(t) = t — c\ In 



1 + — 

ci 



(16) 



Then the two-loop FAPT resummation procedure produces the following answers: 
For the series ()15ap the resummed expression is 



S V [L] = i?x +v [L] 

Cl t 
,Cl +t 

whereas for the series (|15b|) 



+ 



ci+t 

?2+u[L] 



z v dz Ti+u [L + r(tz)~ r(t)] 



1 t 2 z u+i 
dz ■ 



Cl + 1 z 



P(t) 

P 2 +u[L + T{tz)-T{t)] 



(17a) 



p(t) 



[L] 



Bl+vo;vi[L] + 



Cl 



Cl+t 



1-fl 



lp{t) 



Cl t 
(ci + 1) 1 -^ 



dz 



dz 



dB 1+U0 , Ul [L + T(tz)-r(t)] 



(ci +tz) 



"1 



dL 



P(t) 



(ct+tz) 1 



VQ+Vl-l 



B 2+U0 . U1 [L + T(tz)-T(t)\ 



(17b) 



P(t) 



5 



with 5q j1/q+Ui being a Kronecker delta symbol. 



3 Applications to Higgs boson decay 

Here we analyze the Higgs boson decay to a bb pair. For its width we have 

R S (M 2 H ) 



r(H -> bb) = rgK 



3 m 2 



(18) 



with r, 



mt 



3 Gf Mh m^/4-v/2vr, mj, and Mh are the pole mass of the 6-quark and the mass 
of the Higgs boson, respectively, R S (M H ) = m 2 (M H ) R$(M H ) and R s {s) is the i?-ratio for the 
scalar correlator, see for details in [16} 125]. 

In the one-loop FAPT this generates the following non-power series expansion^ 



ft s (1) [L] =3m 



(i) 



n>l 



(19) 



where = ^-53 ^ ^-09 GeV 2 is the RG-invariant of the one-loop m 2 (/i 2 ) evolution, written 
in the following way: m 2 (n 2 ) = mf^ at%° (fi ) with z^o = 2[7o/6o]A r / =5 = 1-04, and 70 is the 
leading-order quark-mass anomalous dimension. This value ih 2 ^ has been obtained using the 
one- loop relation [261 I2Z] between the pole mass of the b quark, m&, and the value of the 
running mass at the scale //* = mj(^), which we call mf,(m 2 ). Here we also extract the value 
of df = 17Cf/4 = 17/3 out of higher perturbative coefficients, so that d% = 3d^/17 and df = 1. 
In the two-loop case we obtain 

ci b a s (/j, 2 ) 



n^ ] [L] = 3m 



m 2 2) <V) 



1 + 



4tt 2 



U,. n 



l n>l 



(20a) 
(20b) 



where the RG-invariant mass rh 2 ^ 
the one-loop case, and v\ 



8.22 ± 0.09 GeV 2 , while the value of vq is the same as in 
2 [(71 bo — 70 61)/ (60 &i)]jv>=5 = 1-86 (71 is the next-to- leading-order 
quark- mass anomalous dimension). Note that we determined both values, mL and mjL, using 



the estimates for mj(mj) derived by Penin and Steinhauser [28] . 

In order to estimate now the importance of higher-order corrections, we need to construct 
some model for the generating function P(t). We use a Lipatov-like model [23] with c = 2.4, 
f3 = —0.52, which generates factorially growing coefficients d^: 



d s 



Ps(t) 



(t/c)+/3 
c(l + /3) 



1 + 

-t/c 



(21a) 
(21b) 



this model gives a very good prediction for d^ with n 



61.1, and 625 in comparison with 7.42, 62.3, and 620. Moreover, it predicts the value di 



2, 3, 4, calculated in QCD PT [25J: 7.50, 

7826 

which is in a good agreement with the PMS [29, 30j prediction d™ s = 7782, obtained in [24J. 

Then, we apply the FAPT resummation technique to estimate how good is the iVth order 
truncation 



in 



(i) 



N 



d;, 



(22a) 



4 Appearance of denominators ir n in association with the coefficients d n is due to d n normalization. 
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in approximating the whole one-loop width r^^ fe -[L] in the range L S [12.4, 13.5] which corre- 
sponds to the range M H £ [100, 172] GeV with Aq£,= 3 = 201 MeV and 2lf lob (m|) = 0.1226. In 
this range, we have L§ < L < Lq = ln(mf/A|) and hence 



p(l);oO rjri 



m 2 {1) 



2CM + f (ai+M, 



vr/3 5 



+ A 6 2l 



(22b) 



pS 
"0 ■ 



with i^ft) defined via Eqs. (jlip and (|2ip . We analyze the accuracy of the truncated FAPT 
expressions (|22ap and compare them with the total sum (|22b|) using relative errors A^[L] = 

2, N = 3, and N = 4 in 



r„ , T \L:N]/T [ L r J ' a !i' T \L]. In Fig. 1 we show these errors for A" 

(1) 



-.(i);oo 



the analyzed range of L G [11, 13.7]. We see that already r^/^ fe -[L;2] gives an accuracy of the 



order of 2.1 



o, whereas ^[L; 3] is of the order of 1%. Looking in the left panel of Fig. 1, we 



.025 


A 2 [L] 




. 02 






0.015 






. 01 


A 3 [L] 
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13 
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Figure 1: Left panel: The relative errors Ajv[L], N = 2,3 and 4, of the truncated FAPT in 
comparison with the exact summation result, Eq. (|22b|) . Right panel: The width T H .j^ as a 
function of the Higgs boson mass Mh in the resummed FAPT. The width of the shaded strip 
is due to the overall uncertainties, induced by the uncertainties of the resummation procedure 
and the pole mass error-bars. Both panels show the results obtained in the one-loop FAPT [24J. 



understand that in order to have an accuracy better than 0.5%, we need to take into account the 
4-th correction. We also verified that the uncertainty due to P(t)-modelling is small, < 0.6%, 
while the m^rn 2 )-induced uncertainty is of the order of 2%. The overall uncertainty then is of 
the order of 3%, see in the right panel of Fig. 1. 

Qualitatively, the same picture is reproduced at the two-loop order [Mj. In the left panel 
of Fig. [21 we discuss the convergence properties of the decay widths, truncated at the order 
N, relative to the resummed two-loop result ™ T (Mj). We see that our conclusions drawn 

H — >bb 

(2) 

from the one- loop analysis remain valid. Indeed, T H -(M#;2) deviates from the resummed 

(2) 

quantity by not more than 2%, whereas F v J -{Mh; 3) reaches an even higher precision level of 
the order of 0.7%. 

In the right panel of Fig. [21 we show the results for the decay width -r(Mff) in the 
resummed two-loop FAPT, varying the mass in the interval rh^) = 8.22 ± 0.13 GeV according 
to the Penin-Steinhauser estimate mfe(m^) = 4.35 ± 0.07 GeV [2H]. Comparing this outcome 
with the one-loop result — upper strip in the same panel of this figure — reveals a 5% reduction 
of the two-loop estimate. This reduction consists of two parts: one part (w +7%) is due to 
the difference — m both approximations, while the other (~ —2%) comes from the 
difference in the values of Rs(Mh) in the one- and the two-loop approximations. 
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Figure 2: Left panel: The two-loop width T H ^ as a function of the Higgs-boson mass Mjj 
in the resummed (black solid line) and the truncated (at the order N) FAPT |24| . Here, the 
short-dashed line corresponds to JV = 1, the dashed one to N = 2, and the long-dashed line to 
N = 3. Right panel: The two-loop width r°° ,r as a function of the Higgs-boson mass Mh in 
the resummed FAPT is shown (lower strip). The upper strip shows the corresponding one- loop 
result. 



4 Application to the vector correlator Adler function 

Now we analyze the Adler function of the vector correlator. We model the generating function 
of the perturbative coefficients d n [31] (note that here d\ = 1, so that d n = d n ) 

D Y [L] = d + jrd n f^BY (23) 

n=l V ^ / 



using a generalized Lipatov-like model: 



xn+1 _ , , 

d n = cn ~ l 5 2_ 1 T (n), (24a) 



5e-^ s -(t/c)e-^ 
P vW " ^T^T) • ( 24b ) 

Our predictions, obtained with this generating function by fitting the two known coefficients 
and ^3, have been included in Table [T] (see [M])- We observe a good agreement between 
our estimate d\ = 27.1 and the value 27.4, calculated recently by Chetyrkin et al. in Ref. 
|31| . Would we use instead a fitting procedure, which would take into account the fourth-order 
coefficient d^ in order to predict d$, we should readjust the model parameters in (f24"j) to the new 
values {c = 3.5526, S = 1.32453}. In order to explore to what extent the exact knowledge of 
the higher-order coefficients d n is important, we employed our model (|24p in [23] with different 
values of the parameters: c = 3.63 and 5 = 1.3231. These values were selected in order to 
reproduce the exact value of the coefficient d^ = 27 A by the value, close to the Naive Non- 
Abelization (NNA) prediction |32| . The difference between the resummed results of the two 
models in the region of Nf = 4 appears to be very small — of the order of 0.2%. We also 
show in Table [1] the d n estimates obtained by using the Improved NNA (INNA) approximation, 
developed in our last paper [Mj by taking into account not only the b^a™ terms, but also the 
bibQ~ 2k ctg ones. This technique produces predictions for higher-order coefficients in terms of 
intervals for the d n values. We see that our model predictions for the d$ coefficient are all inside 
the INNA interval. 
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Table 1: Coefficients d n for the Adler- function series with Nf = 4. The numbers in the square 
brackets denote the lower and the upper limits of the Improved NNA (INNA) estimates. 



PT coefficients 


d x 


d 2 


^3 


di 


d 5 


pQCD results with Ay = 4 [31] 


1 


1.52 


2.59 


27.4 




Model (J23D with c = 3.555, 5 = 1.3245 


1 


1.52 


2.59 


27.1 


2024 


Model (J21J) with c = 3.553, 5 = 1.3245 


1 


1.52 


2.60 


27.3 


2025 


Model dMJ) with c = 3.630, 6 = 1.3231 


1 


1.53 


2.26 


20.7 


2020 


"INNA" prediction of [M] 


1 


1.44 


[3.5,9.6] 


[20.4,48.1] 


[674, 2786] 


Then, we have 












D$ PT (Q 2 ) 




n>l 


At h (Q 2 ), 




(25a) 


D^(Q 2 ;N) 




N J v 

n=l 


At h (Q 2 )- 




(25b) 



The global-APT resummation result for Dy PT (Q 2 ) can be estimated using our resummation 
formulas from [18 1 1231 12"4"] . which we shortly explained in Sect.[2J We show in Fig. [3] the relative 
errors A^(Q 2 ) = 1 - D§ PT (Q 2 ; N)/D§ PT (Q 2 ) evaluated in the range Q 2 G [2,20] GeV 2 for 
four values = 1,2,3,4. We observe from this figure that already -Dy PT (Q 2 ;l) provides 
an accuracy in the vicinity of 1%, whereas D^f 1 \Q 2 ;2) is smaller then 0.1% in the interval 
Q 2 = 1 — 20 GeV 2 . This means that there is no real need to calculate further corrections: at 
the level of an accuracy of the order of 1% it is quite enough to take into account only the 
terms up to d 2 . Again, we verified that this conclusion is stable with respect to the variation 
of the model Py{t) parameters. This conclusion is in some sense surprising: The best order 
of truncation of the FAPT series in the region Q 2 = 2 — 20 GeV 2 is reached by the N 2 LO 
approximation, i.e., by keeping just the c^-term. 
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Figure 3: The relative errors A^(Q 2 ) evaluated for different values of A: N = 1 (short-dashed 
red line), N = 2 (solid blue line), N = 3 (dashed blue line), and N = 4 (long-dashed blue 
line) of the truncated APT given by Eq. (|25bp , in comparison with the exact result of the APT 
resummation procedure. 
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5 Conclusions 



We described here the resummation approach in the global versions of both APT and FAPT 
for the one- and two-loop running of the effective QCD coupling. We argued that this approach 
produces finite answers, provided the generating function P(t) of the perturbative coefficients d n 
is known. We considered two applications of our approach to the physical problems of current 
interest. 

In the case of the Higgs boson decay, H — > bb, the main conclusion is that to achieve a 
truncation accuracy of the order of 1% for the width r^_^(Mff) it is enough to take into 
account up to the third correction — in complete agreement with Kataev&Kim [33]. In the 
case of the vector-current Adler function, T>y(Q 2 ), the truncation accuracy is of the order 0.1%, 
already reached at N 2 LO, i. e., by taking into account up to the second-order correction. 

In our approach the knowledge of higher-order coefficients, d^, is needed in order to construct 
a reliable model for the generating functions P{t). Meanwhile, we showed that the exact values 
of still higher coefficients, d§, are not so important. 
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